nucl-th/9608012 

Quasielastic K + -mic\eus scattering 



A. De Pace 

Istituto Nazionale di Fisica Nucleare, Sezione di Torino, 
via P.Giuria 1, 1-10125 Torino, Italy 



C. Garcia- Recio 

Departamento de Fisica Moderna, Faculdad de Ciencias, 
Universidad de Granada, E-18071 Granada, Spain 



E. Oset 

Departamento de Fisica Teorica and Instituto de Fisica Corpuscular, 
Centro Mixto Universidad de Valencia - Consejo Superior de Investigaciones Cientificas, 

4-6100 Burjassot, Valencia, Spain 
(August 1996) 



Abstract 

Quasielastic -fT + -nucleus scattering data at q = 290, 390 and 480 MeV/c are 
analyzed in a finite nucleus continuum random phase approximation frame- 
work, using a density-dependent particle-hole interaction. The reaction mech- 
anism is consistently treated according to Glauber theory, keeping up to two- 
step inelastic processes. A good description of the data is achieved, also 
providing a useful constraint on the strength of the effective particle-hole 
interaction in the scalar-isoscalar channel at intermediate momentum trans- 
fers. We find no evidence for the increase in the effective number of nucleons 
participating in the reaction which has been reported in the literature. 
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I. INTRODUCTION 



Quasielastic studies are traditionally a good source of information about nuclear and 
nucleon structure. The main tool has been usually represented by electron scattering ex- 
periments, since in this case the elementary probe-nucleon reaction mechanism is regarded 
to be under better control. However, not all the possible nuclear response functions can be 
accessed in this way and, furthermore, the extraction from the data of the response functions 
that enter electron scattering can be model dependent. 

Hadronic probes are an alternative source of quasielastic data. Among them it is par- 
ticularly interesting the case of _fT + -nucleus scattering, since, for kaon laboratory momenta 
below 800 MeV/c, the elementary f^ + -nucleon (K + N) interaction is much weaker than 
other hadron- nucleon interactions. In fact, a weaker elementary interaction allows the pro- 
jectile to penetrate deeper inside the nucleus, thus probing regions of higher density and, 
consequently, being more sensitive to collective phenomena. 

Such an experiment has been performed at BNL, where quasielastic i^ + -nucleus double- 
differential cross sections have been measured for kaons with a laboratory momentum of 705 
MeV/c, using D, C, Ca and Pb as targets 0,0]- Data have been taken for scattering at 24°, 
34° and 43°, corresponding to approximately fixed momentum transfers of 290, 390 and 480 
MeV/c. In the case of C, the data have been taken for all the momenta, whereas for Ca and 
Pb they are available at 290 and 480 MeV/c. 

In Refs. [0,0 the data have been analyzed using various relativistic models to describe 
the nuclear dynamics; on the other hand, the distortion of the incoming kaons has been 
accounted for employing a very simple model, based on the effective number of nucleons 
participating in the reaction, N eS (see Sec. [11 D| ). The "experimental" values found for iV eff 
are ~ 30% higher than the ones calculated in the Glauber theory. This finding is rather 
puzzling, also in view of the fact that multiple scattering theory underestimates the nuclear 
elastic scattering data |||4]], a fact which has been interpreted as a signal of an enhancement 
of the in-medium K + N cross section. As also noted in Ref. [0, the quasielastic cross section 
is, in general, proportional to the elementary single differential K + N cross section and to the 
effective number of participating nucleons, the latter depending, in turn, on the total K + N 
cross section (see Sec. |II D| for details). At the energy of the experiment, K + N scattering is 
nearly isotropic, so that the differential and total cross sections are proportional; however, an 
increase in the K + N cross section reduces N e g (and viceversa), leaving a quasielastic cross 
section little dependent on the input K + N amplitudes. Thus, the quasifree experiment on 
the one side does not shed light on the elastic scattering data and, on the other, it seems to 
pose another puzzling problem. 

In this paper we would like to reanalyze the i^ + -nucleus quasielastic data, using a more 
realistic model for the reaction mechanism, — already applied to (p, p') and (p, n) quasifree 
scattering - - based on a consistent implementation of Glauber theory for quasifree 

scattering. The nuclear dynamics has been treated through a finite nucleus continuum 
random phase approximation (RPA) calculation, also accounting for the coupling of the 
particle- hole (ph) states to higher order configurations (spreading width). The dynamics is 
non-relativistic, but the correct relativistic kinematics has been used. 

We will show that the puzzling experimental outcome for N e g is very likely an artifact 
of the naive fitting procedure employed in the analysis of the data and that a realistic (and 
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without free parameters) theoretical framework has no difficulties in giving a rather accurate 
description of the experimental cross sections. 

In Sec. [TI| we introduce the theoretical methods employed in the calculations: Here, one 
can find a brief description of the quasielastic nuclear responses in finite nuclei; of the RPA 
using density-dependent interactions; of the specific model for the ph effective potential that 
we have used; and of the response functions to hadronic probes up to two-step processes in 
the Glauber framework. In Sec. ID, we display the results of our calculations and in the last 
Section we discuss their implications for the interpretation of the data. 

II. THEORETICAL METHODS 

A. Quasielastic nuclear response 

The nuclear response function to an external probe, transferring momentum q and energy 
u, is proportional to the imaginary part of the polarization propagator 0: 

R a {q, to) = --lmU a (q, q, u) . (2.1) 

7T 

The latter reads 

n a {q,q';u) = T,^o\d a (q)\^ n )^ n \dl(q')\^ ) 

1 

X 

huj - [E n - E ) + 17] 

1 

fko + (E n — E ) — ir\ 

where {l^n)} is a complete set of nuclear eigenstates of energy E n , O a (q) the second quan- 
tized expression of the vertex operator and a labels the spin-isospin channel. In the following, 
for the purpose of illustration, we shall deal explicitly only with the scalar-isoscalar channel, 
- which is, by the way, the dominant one in K + N scattering at the energies of concern to 
us, — where 0(q, r) = exp(iq ■ r): The case of the spin modes, which is slightly complicated 
by the spin algebra, is treated in detail in Ref. |5|. 

The angular part of II(qr, q f ; u) (we now drop, for simplicity, the channel label a) can be 
handled through a multipole decomposition that reads 

1%, q', uj) = Y, q'; u)Y$ M {q)Y JM {q') , (2.3) 

JM 

so that 

R(q, u) = -^lmJ2(2J + l)Tlj(q, q, u) . (2.4) 

In a mean field (shell model) framework, one has for the particle-hole (ph) polarization 
propagator: 



(2.2) 
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n° J (g,g'^) = EQi J0 (?) 

ph 



Qi J0 V) 



u - (e p -e h )+ir) to + (e p - e h ) - irj 
where p (h) labels a complete set of single particle (hole) quantum numbers and 

QifiQ) = (jph; J\l*; J)(-iY +1 (-iY h 2[M^ P + i)(24 + i)] 1/2 x^(<?) ( 

In (|2.6|) , (jpjh] J\£cr; J) is the standard LS — jj recoupling coefficient and 

Ze, P h(q) = / drr 2 j £ (qr)R p (r)R h (r) , 

J 



(2.5) 



i P 4 ^ 





(2.6) 



(2.7) 



Rp{h)( r ) being the radial particle (hole) wave function and e p ^ the associated eigenvalues. 
They are obtained by solving the Schroedinger equation with the Woods-Saxon potential 



W(r) 



1 _|_ e (r-R)/a 



+ 



he 

mlc 2 



1 2 



,{r-R)/a 



ar [1 + e (r-R)/aV 



(2.8) 



(neglecting, for simplicity, the Coulomb term) , where is the pion mass and the following 
set of parameters has been employed: 



W = -54.8 MeV, W so = -lOMeV 
R = 1.27 A 1 / 2 fin, a = 0.67 fm . 



(2.9) 



Note that the sum appearing in Eq. fl2.5| ) should be understood as a sum over the discrete 
part of the spectrum and an integration over the continuum one. Indeed, contrary to the 
widespread procedure of calculating the polarization propagator in coordinate space and 
then Fourier transform to momentum space, we have directly evaluated Eq. ( |2.5| ) in the 
latter. Besides being fast and reliable, this procedure allows also for an important extension 
to the shell model polarization propagator, namely the inclusion of the spreading width of 
the ph states. 

This can be accomplished by adding in Eq. (|2.5| ) a complex ph self-energy, i. e. through 
the following substitution: 



(e P - eft) -> (e p - e h ) - T, ph (u) . 



(2.10) 



Although Spft could in principle be calculated, we shall actually employ a phenomenological 
parameterization, writing 



£ p /i(u;) 



A ph (u) + i- 



r p h(o;) 



(2.11) 



with 



T ph (uj) = [%(Huj + e h ) + 7ft(e p - fw)](l + C TS ) 
\ h (uo) = [A p (huj + e h ) + A h (e p - hu)](l + C TS ) 



(2-12) 
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The arguments of the functions ^pih) and A p (m have been inferred from the analysis of 
the second order particle and hole self-energy contributions to E p ^. The coefficients Cts 
represent the corrections, in each isospin T-spin S channel, due to the ph interference 
diagrams: These have been estimated || to contribute around 5% in S = T = 1 and 
~ 20 -T- 25% in S — 0, T — 1 and S — 1, T = 0; given the smooth dependence of the 
quasielastic response on the ph self-energy and the small contribution of these modes to the 
A^-nucleus cross section, these corrections are not important and we shall set for simplicity 
Cts — in these channels. Concerning the channel S — T — 0, the estimate of Ref. || gives 
Coo = ~ 1) resulting in a complete cancellation of the ph spreading width. This outcome is 
valid in nuclear matter and in the limit of very large angular momenta, which is presumably 
good in the quasielastic peak (QEP) region. However, at low transferred momenta the region 
of small energy transfers of the quasielastic response is dominated by resonances associated 
to specific (low) angular momenta. In that case it is a better approximation to retain the 
full width and, accordingly, we shall choose C o = — 1 + exp[— (u/u ) 2 } with u Q = 30 MeV: 
Then, resonances in the 10 -j- 20 MeV range get a width of a few MeV, whereas for u ^ 30 
MeV there is essentially no spreading. 

Finally, the J p (h) are chosen according to the parametrization of Ref. ||, namely 



7 P (e) 



7ft (e) = 7 P (- 




(2.13) 



symmetrical with respect to the Fermi energy (ep = 0), which gives a reasonable fit of the 
particle widths for medium-heavy nuclei, using a = 10.75 MeV, eo = 18 MeV and e\ = 110 
MeV 0. The corresponding real parts are obtained through a once-subtracted dispersion 
relation 0. 

A final remark, concerning relativistic kinematics effects on the response functions, is 
in order. We shall be concerned with momentum transfers up to 500 MeV/c, where purely 
kinematical relativistic effects are starting to be sizable. At q — 500 MeV/c the relativistic 
position of the QEP is ~ 8 MeV below the non-relativistic one, the effect being even larger 



for energies on the right hand side of the QEP. In Ref. | 10| , it had been shown that the non 



relativistic Fermi gas can be mapped into the relativistic one through the simple prescription 



UJ 



+ (2.14) 



The validity of ( 2.14p in a finite nucleus shell-model calculation has been checked in Ref. [JTT 



The main effect of properly accounting for the relativistic dispersion relation turns out to 



be a moderate shrink of the response functions at the right of the QEP (see Sec. |T|). 



B. RPA with density-dependent interactions 

The polarization propagator introduced in the previous subsection is used in this paper 
as input in a continuum RPA calculation of the response functions. 
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For a central density-independent interaction one has to solve, for each multipole J, 
the following integral equation (see Ref. [12] for the complications introduced by the tensor 
interaction): 

nf A (g,g» = n J (g,g» + ^-p j / o dkk 2 Tl°j(q,k;u)V(k)Tl^ PA (k,q';uj) , (2.15) 

V(k) being the Fourier transform of the two-body potential. However, realistic effective 
nuclear interactions are in general density-dependent, especially in the scalar channels. 

For a general interaction potential V(r, r'), one can introduce the double Fourier trans- 
form 

V{k,k') = J drdr'e- lk r e ik '- r 'V{r,r') 

= £ Vj(k, k') Yjm (k) YJ M (k) . (2.16) 
j 

Then, instead of Eq. (|2~T5p , one has to solve the following RPA equations: 

I poo roo 

nf A (q, q'; U ) = H°j(q, q'; u) + -—-^ J q dk k 2 ^ dk' k U j(q, k; co)V K (k, k')Tiy A (k>, q'; U ) 
= U j(q, q'; u) + — ^ dk' k' 2 Kj(q, k'\ w)nf \k', q'; u) , (2.17) 
having defined the kernel 

Kj(q,k';u) = — j Q dkk 2 U j(q,k;u)Vj(k,k') . (2.18) 

Clearly, in the case V(r,r') = V(r-r') one has Vj(k, k') = {2TifV{k)5{k~k')/k 2 : The ker- 
nel ( ^.18[ ) is then reduced to Kj(q, k'\ u>) = Hj(q, k'; u)V(k') and one gets back to Eq. ( |2.15| ). 

From Eq. ( |2 . 1 7 ) it is apparent that the solution of the RPA equations with density- 
dependent forces does not pose any additional technical problem, apart from the input 
kernel, whose calculation is more involved. 

As discussed in detail in the next subsection, in the calculations of Sec. |II|we have been 
using a parameterization of the effective nuclear interaction linear in the density. In general, 
in any (non-tensor) channel one has 



V( ri , r 2 ) = V cx ( ri - r 2 ) + V"(n - r 2 ) p 



ri + rr 



= V ex (r) + V p (r) p(R) , (2.19) 

with r = ri — r 2 , R = (r 1 + r 2 )/2 and p(R) = p(R)/p(0), p(R) representing the nuclear 
density, here approximated by the standard Fermi distribution. 
In momentum space (|2.19| ) reads 

V(k, k') = V ex (k)(2ir) 3 5(k - k') + V P ( 1 1±K] p{k _ k') (2.20) 
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and its Jth component in the angular momentum expansion (|2.16 ) turns out to be 



Vj(k, k') = V ex (k)^6(k - k') + Vj p (k, k') (2.21) 



where 

V^)^J d{k .k)v^^) Kk - M .k)^ (2,2, 
Pj representing the ordinary Legendre polinomials. 

C. Effective ph interaction 

Two strategies are possible in order to determine the effective ph interaction in the 
nuclear medium: One can either directly fix an effective potential by fitting some phe- 
nomenological properties or start with a bare nucleon-nucleon interaction and calculate the 
related G-matrix. Parameterizations of the ph interaction based upon the first procedure 
are generally only available at very low momentum transfers (in terms of Migdal-Landau pa- 
rameters); since we are probing relatively high momenta, we have resorted to use a G- matrix 



and we have chosen the calculation of Ref. ||13|| , which, in our view, has the following ap 



pealing features: It is based upon a realistic boson exchange potential; (nonlocal) exchange 
contributions are included in the effective interaction, which is conveniently parameterized 
in terms of Yukawa functions; a parameterization of the density dependence is also provided. 

The potential is given using the standard representation in spin and isospin (no spin-orbit 
contribution will be considered in the following): 

V{k f , k F )=F + F't x ■ t 2 + Go x ■ a 2 + G'ct 1 ■ cr 2 n • r 2 

+TS 12 (q) + T'S 12 (q)r 1 ■ r 2 + HS 12 (Q) + H'S 12 {Q)r l ■ r 2 , (2.23) 

where q = ki — kf, Q = ki + kf (ki, kf being the relative momenta in the initial and final 
state, respectively) and the coefficients are density and momentum dependent. 

However, before utilizing the nuclear matter interaction of Ref. |13] in a finite nucleus 
calculation of quasielastic responses, a few issues have to be addressed. 

a) The density dependence of the G-matrix is given in terms of density dependent cou- 
pling constants, which is not very useful for our purposes. Furthermore, the parameterization 
is fitted for 0.95fm _1 < kp < 1.36fm _1 , kp being the Fermi momentum, which spans a range 
of densities down to roughly 1/3 of the central density: Extrapolation of their parameteriza- 
tion to lower densities gives unreasonable results. Thus, we have chosen to employ the linear 
p dependence of Eq. (|2.19| ) [or ( |2.20| )1, which is known to provide a reasonable parameteriza- 



tion (see, e. g., Ref. [13]). One can see in Fig. |T] a comparison of the two parameterizations 
for the kp dependence of the effective interaction. It should be noted that most of the 
contribution to the quasielastic responses comes from densities where the two descriptions 
differ by a few per cent. 

b) In order to obtain, at a fixed density, a local interaction, one can use the relation 



between q, Q and ki, i. e., Q = J4:kf — q 2 , substituting ki with a suitably chosen average 
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FIG. 1. Effective interaction in the non-tensor channels as a function of hp at q = 0; linear 
(solid) and from Ref. |l3| (dashed) density dependence. 
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FIG. 2. (a) Direct and exchange ph matrix elements; (b) direct pp matrix element. 
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value, (ki); then, the only independent momentum is q. The authors of Ref. [|13 



were 



interested in a potential for nuclear structure calculations: Hence, they put the momenta 
of the two nucleons in the initial state on the Fermi surface and averaged over the relative 
angle, getting (ki) ~ 0.7 hp- Clearly, in this case one has the constraint < q < lAkp- 
On the other hand, we are interested in the ph interaction in the quasielastic region: One 
nucleon in the initial state is below the Fermi sea, while the other can be well above it. 
A look at Fig. |2] shows that fc, is defined, in terms of the particle and hole momenta, as 
ki = (p — h')/2 = (h — h' + q)/2. Thus, at fixed q one should average ki over h and 



h', getting (ki) « J6kp/5 + q 2 /2. Now ki is growing with q, so that there are no longer 



constraints on q, and the exchange momentum turns out to be constant, Q = y6/5kp. One 
can see in Fig. |3] the resulting interaction in the non tensor channels. 

c) In the case of the tensor channels things are simpler, since in the parameterization of 
Ref. there is no explicit density dependence (Fig. The coefficients of the exchange 
tensor operator, H and H', display a very mild density dependence, induced by Q, which 
is completely negligible. The only catch concerns the treatment of ^12 (Q): Assuming that 
q and Q are orthogonal, with some algebra one can show that Su(Q) = — Su(q)/2- Note, 
however, that 80 -r- 90% of the quasielastic cross section for K + scattering is due to the 
scalar-isoscalar channel. 



D. Response functions with hadronic probes 

The formalism introduced above is not enough when one is dealing with strongly inter- 
acting probes, in which case a framework for the reaction mechanism must be provided. For 
this purpose, we have chosen the Glauber approach [fL5H , including up to two-step inelastic 
processes. 

Simple treatments within the Glauber theory usually amount to including the effects 
due to rescattering in an effective number of nucleons participating in the reaction, thus 
effectively renormalizing the response functions defined above. In the case of K + quasielastic 
scattering, this approach has been followed in Refs. P|J2]JTB[]. 

However, a consistent treatment within Glauber theory leads one to the definition of 
surface response functions [j|]. A detailed derivation, in particular for the spin-isospin chan- 
nel, can be found in Ref. ||. Here, we briefly sketch the relevant points, again using the 
scalar-isoscalar channel as example. 

A surface polarization propagator can be obtained from (|2.2j ) by substituting the vertex 
operator O a (q,r), which describes the probe-nucleon coupling, with 

0™%, r) 1 



Wf a (q) 
x J dbd\e-^ Tib)/2 e l( - q ' x > b f a (X) 
xO a {X,r) , (2.24) 

b and A being bidimensional vectors in the plane orthogonal to the direction of motion of the 
projectile, f a the elementary probe-nucleon amplitudes, a tot the effective total probe-nucleon 
cross section and 
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T(b) = / + °° dzp{r = Vb 2 + z 2 ) , 

J — CO 



v 2 -25) 



p{r) being the nuclear density. With these definitions, the effective number of participating 
nucleons is given as 



N cS = / dbT{b)e 



-<T to tT(6) 



^2.26) 



cf to t depends on the energy and we shall take the average over the total cross sections at 
the incoming and outgoing projectile energies: For K + scattering at k — 705 MeV/c, the 
variation of a tot is rather weak, going roughly from 14 mb to 13 mb for energy losses up to 
250 MeV. 

One then finds, for the J-th multipole, 



n^)(g,g;w) = n J(a) (g,g 
1 



having set 



dXX dX'X'Re[r a (X)f a (X')Gf(X,X'-q)}n J{a) (X,X';u;) 

\ja{q)\ JO Jo 



1 



\Uq)\ 2 J° 



dX A Re[f* a (q) f a (X) Hf\\; q)]U J(a) (q, A; u) , 



Gf(X, A'; q) = Yl CJrn9* m (\ q) 

m 

Hj ] (X;q) = ^2c Jm g m (X,q) , 



(2.27) 

(2.28a) 
(2.28b) 



where 



and 



g m (X, q)= dbb\l- e-^ T ^ 2 ] J m (Xb) J m (qb) 
Jo 1 J 



(J-m- l)!!(J + m- 1)!! 
(J + m)\\(J - m)\\ 



' J+m 



0, J + m odd 

1, J + m even 



(2.29) 



(2.30) 



With the previous definitions, the one-step surface response functions are given as 
i?™ rf (g, a;) = — Im X)j(2 J + l)nj urf (g, q; u)/47i 2 and the double differential cross section turns 
out to be 



d 2 a 



dVtde' 



= £l/«(?)l 2 K ur W), 

l-stcp a 



(2.31) 



to be compared with the effective number approximation, which reads 
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d 2 a 



dVtde' 



N cS J2\fa(q)\ 2 R a (q,uj) . 



(2.32) 



l-stcp 



Multiple scattering inelastic contributions can, in principle, be incorporated along the same 
lines. Since in the kinematic regions of interest they turn out to be much smaller than 
the one-step terms, it is not worth going through a very complex formalism and one can 
safely stick to the effective number approximation. In this case, the two-step contribution is 
proportional to the convolution of two one-step response functions [HJ . By summing over all 
the possible channels and using free one-step response functions, one can define the two-step 
response as 



1 



q,uj) 



dq' / dJ 

h 



k 2 EsT\fsAq)\ 2 J Jo 
fe |/st(?0| 2 # (O) (?V) E \fsA\q ~ q'\)\ 2 R ( °\\q -<J\,u>- «0 



ST 



ST 



+2Y / \fsi(q')\ 2 R { °\q',u;>) £ \f sl (\q - q'\)\ 2 R i0) (\q -q'\,u- J) 

s s 



[2M) 



where k is the momentum of the projectile, fsr the KN amplitude in the S, T channel and 



U> ^ I dbT 2 (b)t "■ 



[2.U) 



is proportional to the effective number of pairs participating in the double scattering. 

The full quasielastic fC + -nucleus cross section discussed in Sec. |II| is given by the sum 
of one- and two-step terms: 



d 2 a 



dQde' 



7 = Ei^(^)i 2 ^r f (?^) + J R (2) ('?^) 



(2.35) 



We conclude this Section with a few remarks about the elementary K + N amplitudes 
employed in the calculations. The amplitude for elastic K + N scattering can be written as 
yO) _|_ f^r 3 , with /M = g(v -f icr ■ hh^\ n being a unit vector normal to the scattering 
plane. In Fig. |5| one can see the amplitudes as a function of the momentum transfer for 
k = 705 MeV/c [|l^,[n|: The dominant channel is clearly the scalar-isoscalar one. 

Some care should be taken in the choice of the reference frame where the K + N amplitudes 
are evaluated. In order to be able to factorize the two-body amplitude out of the nuclear 
response functions, an "optimal" choice of the reference frame has to be done |]20|j2l[] : The 
optimal momentum of the struck nucleon turns out to be 



(2.36) 




such as to reduce to zero at the QEP and to the Breit frame value at uj = 0. 
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FIG. 5. Center-of-mass K + N squared amplitudes corresponding to a K + laboratory momen- 
tum k = 705 MeV/c. 



The most relevant consequence of this choice of frame is that the elementary projectile- 
nucleon amplitudes should be evaluated at an effective laboratory kinetic energy defined 
by 



neflf 



E k E opt - k ■ p opt - m p m N 



m N 



(2.37) 



m p and k being the projectile mass and momentum, respectively. This introduces a depen- 
dence on u of the amplitudes entering in the evaluation of the quasielastic cross section, 
which can be important if the former are strongly energy dependent: At the kinematics of 
relevance to us, the effective K + momentum varies roughly in the range 500 4- 800 MeV/c, 
where the amplitudes change sufficiently slowly to make the effects not at all dramatic. 



III. RESULTS 

Let us start by briefly discussing the relativistic kinematical effects, introduced at the 
end of Sec. [11 A| . In Fig. ^| one can see the RPA scalar-isoscalar response of 12 C to K + 
probes in the non-relativistic case compared to the response function where the relativistic 
dispersion relation has been accounted for. As expected, the effect is negligible at the lowest 
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FIG. 6. Scalar-isoscalar RPA response function for K + quasielastic scattering: Non-relativistic 
(dotted) and with the kinematical relativistic corrections discussed in the text (solid). 



momentum transfer (q = 290 MeV/c), whereas at higher momenta (q = 480 MeV/c) it 
produces a moderate shrink of the response at the right of the QEP, compensated by an 
enhancement (~ 10%) at the peak position. 

Next, we would like to discuss the uncertainties connected to the choice of the effective 
ph interaction for quasielastic calculations. The knowledge of the latter suffers, in general, 
from many shortcomings: When it is theoretically calculated, — as the one we employ here 
(see Sec. |II C\) , — sources of uncertainties arise because of the specific many-body scheme 
which has been adopted, from relativistic effects (see Ref. [^2j for an effective interaction 
along the lines of Ref. but in a relativistic context) and from the approximations made 
in the actual calculations or in extrapolating to the ph kinematical regime (see Sec. |11 C| ); 
when the effective interaction is just fitted to phenomenological properties, it can suffer from 
ambiguities in the parameterization (many sets of parameters reproducing the same body of 
data) and from the limited range of momenta which is covered (see, for instance, Ref. |14j 



for a phenomenological density-dependent interaction expressed in terms of Migdal-Landau 
parameters at q ~ 0). 

On the other hand, one can of course reverse the argument and use the quasielastic studies 
to gain insight into the effective ph interaction at momentum transfers of a few hundreds 
MeV. One obvious difficulty in this case is related to the fact that it is not always possible 
to disentangle the various spin-isospin channels: The only reaction for which this has been 
achieved is the (p, ft) charge-exchange one, where the separated isovector spin-longitudinal 
and spin-transverse responses have been extracted P3| . However, the strong interaction 
of protons with nuclei constrains that reaction to the low density peripheral region of the 
nucleus, making it little sensitive to RPA effects; moreover, while the spin-longitudinal 
channel can be well described, the transverse one shows much more strength than expected 
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FIG. 7. K + - 12 C quasielastic response functions R S t (S,T = 0, 1) at q = 480 MeV/c. The 
dotted line represents the uncorrected case; the solid line the RPA responses corresponding to the 
interaction of Sec. II C ; the dashed lines the RPA responses with the same interaction scaled by 
±50%. 



||. In this regard, the i^ + -nucleus reaction is much more promising for two reasons: First of 
all, the small K + N cross section makes the kaon enter inside the nucleus much more deeply 
than protons or pions, in turn implying stronger collective effects; secondly, as seen in Fig. |5], 
the scalar-isoscalar channel is largely dominant, making this reaction a "quasipure" probe 
of the S = T = mode. 

Whatever be the attitude towards this problem, we believe it is useful to gain some feeling 
about the sensitivity of the collective effects to the input effective potential. For this purpose, 
we display in Fig. [7|, for the four spin-isospin channels, the surface RPA response functions 
at q — 480 MeV/c, calculated with the ph interaction discussed in Sec. [11 C| , comparing them 
to the RPA responses calculated with the same interaction scaled by ±50%. This amount 
for the rescaling is not an estimate of the theoretical or phenomenological uncertainty of the 
ph interaction; however, concerning the scalar-isoscalar channel, there are reasons to believe 
that the G-matrix estimate of Ref. gives too much attraction: As discussed by the 



authors, the extracted /o Landau parameter makes nuclear matter unstable, a shortcoming 
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that can be cured by higher order contributions^. Thus, while the case with the interaction 
increased by 50% is taken mainly for illustrative purposes, the case with a weaker interaction 
is somewhat more realistic. 

By inspecting Fig. [7] one can note the following: 

i) the dependence of the responses on the potential is somehow desensitized by the RPA 
chain: For the S = T = mode the variation of the response function at the left of the QEP 
is typically about 20 -f- 30%, apart from the very low energy region (u ^ 30 MeV), where it 
grows up to ~ 50%; for the other modes, still at the left of the QEP, the variation is always 
much smaller, say < 5 10%; 

ii) at the right of the QEP the effect of correlations is much smaller, being limited to a 
few per cent for the S = T = mode and to practically nothing for the other modes; 

iii) the relatively weak sensitivity of the isovector and spin channels to variations of 
the input potential makes the small (10 4- 20%) contamination from these modes in the 
if + -nucleus cross section rather stable with respect to uncertainties due to correlations, 
enhancing the argument in favour of this reaction as a probe of the scalar-isoscalar channel. 

The results if Fig. ^ are given for q = 480 MeV/c: Similar considerations apply also to 
the lowest momenta. 

In Fig. |8] one can see the data for the quasielastic K + -C cross section [lljfl compared 
to our calculations. In the theoretical cross sections, besides the case with free response 
functions (dotted), there are included the nuclear responses corresponding to the interaction 
of Sec. |II Q (solid) and to a reduction of 50% of the same interaction (dashed); the curves 
include the two-step contribution, which is also shown separately (dot-dashed). 

From inspection of the figure, one sees that the strength of the cross sections on the 
right hand side of the QEP is well reproduced by both the correlated and the free responses 
(apart from the very high energy tail) at all the momentum transfers; however, at the lowest 
momenta the data show a clear distortion of the typical quasielastic shape, which is not 
compatible with the uncorrelated cross sections. 

The experimental energy resolution is not sufficient to cut out the elastic scattering 
contribution; moreover, low-lying discrete nuclear excitations are much affected by the detail 
of the nuclear model (such as shell model parameters and spreading width). If, for these 
reasons, one excludes the very low energy tail (say, u ^ 15 MeV) it appears that the model 
with a weaker ph interaction gives a remarkably good description of the data at q = 290 
and 390 MeV/c. 

At the highest momentum, the elastic contamination and the discrete excitations have 
been washed out; the interaction in both the correlated models is sufficiently weak to cause 
little distortion of the quasifree shape, but the low energy part of the spectrum seems again 
to favour the model with the weaker ph interaction. 

It is worth noticing the smallness of the two-step contribution, which is at most a few 
per cent of the one-step term at q = 480 MeV/c: This can be contrasted, for instance, to 



Also relativistic effects give rise to a weaker attraction: In Ref. |22| the strength in the S = T = 
channel has been shown to be reduced by nearly a factor 1/2 at q = 0; however, in the range of 
momenta of interest to us, the relativistic and non-relativistic G-matrices appear to be comparable. 
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FIG. 8. Cross section for K + - 12 C quasielastic scattering: Free response (dotted), RPA with 
the interaction of Sec. [I C| (solid); RPA with the same interaction reduced by 50% (dashed). In all 
the cases, the two-step contribution has been added and it is also shown separately (dot-dashed). 
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the isovector spin-transverse response to 500 MeV protons, where one finds, at q = 500 
MeV/c, a sizeable 20 -r- 30% contribution from two-step processes 0. Protons of 500 MeV 
of kinetic energy have a squared momentum which is roughly twice the one of 700 MeV/c 
kaons; furthermore, the two-step factor V 2 of Eq. ( |2.33| ) and ( |2.34j ) is larger by roughly a 
factor two for the kaon reaction, — because of the smaller probe-nucleon total cross section 
c^tot, — but this is compensated by the factor two arising from the two possible orderings 
of the charge-exchange reaction (see formula (8) of Ref. 0). Hence, apart from a factor 2 
in the coefficient, the relative weight of the two-step term in the two reactions has to be 
driven by the elementary amplitudes entering (|2.33j ) (the response functions being in all the 
cases the free ones), that is by the factor |/oo(|<7 — 9 / |)| 2 |/oo(q , ')| 2 /I/oo(q , )| 2 m (P-33|) (since the 
scalar-isoscalar channel is dominant) and by the factor \foo N (\q — q'\)\ 2 \fT N ( c l')\ 2 /\fT N ( c l)\ 2 
in the analogous expression for the (p, n) reaction, q being the external momentum and q' 
the integration variable, whereas f^ N and f^ N are the isovector spin-transverse and the 
scalar-isoscalar NN amplitudes, respectively. 

The different size of the two-step contributions is indeed determined by the different 
strength and the different momentum behaviour of the above amplitudes: In the kaon case, 
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the scalar-isoscalar amplitude has little momentum dependence, so that \foo(q')\ 2 /\foo(q)\ 2 is 
roughly between one and two, while the remaining squared amplitude is on average ~ 1.5 mb 
(in the laboratory); on the other hand, for the (p, n) reaction at q = 500 MeV/c, the ratio 
\f<r N (q') \ 2 l\fr N {l)\ 2 ls typically in the range 4 4- 8, whereas the remaining amplitude gives 
contributions mainly in the range from 10 to 22 mb. In qualitative words, we can say that 
because the NN cross section is relatively forward peaked, — unlike the K + N one, — a 
double scattering in which the momentum transfer q is shared by the two nucleons is favored 
in the NN case with respect to the K + N one. 

Finally, we display also the available data for Ca at q = 290 and 480 MeV/c (Fig. 
To these figures one can apply the same considerations made for the carbon data, although 
the preference for the weaker ph effective potential is even more clear now, since in heavier 
nuclei collective effects tend to be stronger. Indeed, the sensitivity of the response functions 
to variations in the input potential is a little more pronounced in 40 Ca than in 12 C, being 
about 25^35% for a 50% change of the potential in the S = T = channel, to be compared 
to the 20 4 30% sensitivity found in 12 C (see the discussion of Fig. [7]). Also the two-step 
contribution is, as expected, larger in 40 Ca than in 12 C. 

IV. CONCLUSIONS 

From the discussion in the previous Section, it appears that the model presented here 
gives a good description of the quasielastic K + data, the main uncertainty being related to 
the strength of the effective ph interaction in the scalar-isoscalar channel. A good description 
of the data was achieved by using the interaction of Ref. [O] quenched by 50%. It was shown 
that the use of this interaction with the full strength led to results in clear contradiction 
with the data. Indeed, there are theoretical indications that the strength of this interaction 
in the scalar-isoscalar channel should be reduced. On the other hand, from Figs. § and |9| 
one can see that the K + data can be used to constrain the strength of the ph potential for 
the S = T = mode. 

There are, however, two issues that need to be commented upon, namely the reported 
high experimental values for the effective number of nucleons A e g |I|,0] (see the Introduction) 
and the comparison with the available calculations for this reaction that employ relativistic 
dynamical models HI],|2|,^6| . 

Concerning the first point, it might appear curious that the "experimental" value for 
N e g quoted in Refs. [0,0 be ~ 30% higher than the one calculated in the Glauber model, 
since we have seen that the latter gives a nice estimate of the quasielastic cross sections. 
A source of error may of course be given by the use of the effective number approximation 
( [2.32j ); however, we believe that the main reason for the overestimate of N e s lies in the way 
followed to extract the effective number from the data |fj|. Indeed, N e g has been obtained 
by integrating the quasielastic data over the whole range of transferred energies: In order to 
do this, the data have been fitted with Gaussian distributions and, although a subtraction 
of the elastic scattering and low energy excited levels contributions has been attempted, no 
subtraction for any high energy background has been included. From Figs. [| and |9| it appears 
that the very high energy tail of the cross sections is underestimated by the calculations, both 
for the free and the correlated models (remember, from the discussion of Sec. |II|, that while 
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RPA correlations sizably affect the low energy side of the QEP, they have little influence 
on the high energy part). The data actually seem to show some structure beyond the ph 
response region: In the case of quasifree electron scattering, the strength in this region is 
attributed to effects beyond RPA (such as meson-exchange currents or 2p-2h excitations). 
Note also that a simple, "straight line", estimate of the background brings the value of N e Q 



in close agreement with the Glauber estimate p4 | 



Meson-exchange currents effects tied to the interaction of the kaon with the nuclear pion 
cloud have been evaluated in Refs. [p5| , p6| and found to lead to small corrections to the 



i^ + -nucleus total cross section. The consideration of contact terms using chiral Lagrangians 



makes the corrections even smaller |27] , so that this source of corrections cannot account for 
the strength at large values of to. On the other hand, the tail of the A excitation reaches this 
region, since the A acquires a finite width, even below the pion production threshold, due 
to the AiV — > NN reaction in the nucleus. Hence, the channel of kaon induced A excitation 
in nuclei should become a target of both theoretical and experimental investigation in order 
to further clarify these issues. 

Concerning the calculations with relativistic models of Refs. p|j2]Jl6[1, it is obviously 
difficult to make a comparison with our results, since in those papers the RPA response 
functions (calculated in a variety of models) have been multiplied by the "experimental" 
iV e ff, so that their strength is just fitted to the data (of course, the use of the calculated 
value for iV eff would result in a general underestimate of the cross sections). However, it 
should be noticed that the relativistic RPA isoscalar response is quenched, in contrast to 
our non-relativistic RPA response, which is mainly shifted to the left and enhanced at very 
low energies. Then, it is clear that this quenching has to be compensated by a higher N e g. 
Also to be noticed is the fact that the relativistic RPA calculations are able to describe 
(albeit through a fit of the total cross section) only the data at high momentum transfers. 
At q = 290 MeV/c, the large enhancement in the cross section visible at moderate energy 
transfers (15 < c<j < 40) is clearly not predicted, even in a finite nucleus calculation JIJ, in 
contrast to our non-relativistic results. 

Of course, before drawing more firm conclusions one should test the model also against 
the (e, e') data and calculations in this direction are in progress. In this connection, it is 



interesting to note that a recent reanalysis of the (e, e') world data [28 seems to rule out 
the long-standing problem of the missing strength in the charge channel, which had been 
interpreted in terms of a quenching of the charge response. 
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